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Abstract 

Different perturbation theory treatments of the Ginzburg-Landau phase 
transition model are discussed. This includes a criticism of the perturbative 
renormalization group (RG) approach and a proposal of a novel method provid- 
ing critical exponents consistent with the known exact solutions in two dimen- 
sions. The usual perturbation theory is reorganized by appropriate grouping of 
Feynman diagrams of (f 4 model with 0(n) symmetry. As a result, equations for 
calculation of the two-point correlation function are obtained which allow to 
predict possible exact values of critical exponents in two and three dimensions 
by proving relevant scaling properties of the asymptotic solution at (and near) 
the criticality. The new values of critical exponents are discussed and compared 
to the results of numerical simulations and experiments. 

Keywords: Ginzburg-Landau model, Feynman diagrams, renormalization group, 
critical exponents, quenched randomness. 

1 Introduction 

Phase transitions and critical phenomena is one of the most widely investigated 
topics in modern physics. Nevertheless, a limited number of exact and rigorous re- 
sults is available p|. Our purpose is to give a critical analysis of the conventional 
approach in calculation of critical exponents based on the perturbative renormal- 
ization group (RG) theory j|, [|, ||] and to propose a new method which provides 
results consistent with the known exact solutions. The usual RG theory is based on 
several assumptions which could seem to be plausible since the predicted values of 
critical exponents are well confirmed by some numerical results, particularly, by the 
estimations of the high-temperature series expansion [pL Q. The basic hypothesis 
of RG theory is the existence of a certain fixed point for the RG transformation. 
The usual RG theory treatment of the Ginzburg-Landau model is based on the di- 
agrammatic perturbation theory (Feynman diagrams). Straightforward application 
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of the perturbation theory near criticality appears to be problematic in the case of 
the spatial dimensionality d < 4 because of the infrared (i. e., small wave vector k) 
divergence of the expansion terms. Wilson and Fisher |2], |?]] have proposed a way to 
overcome this difficulty by expanding the Feynman diagrams of renormalized per- 
turbation theory in double series of e = 4 — d and In A: (as regards just the critical 
surface). From this the famous e-expansion of critical exponents has originated. 
The first results have been obtained in [0, ^, [J. Nowadays, an explicit e-expansion 
is available up to the fifth order |l(], 11, 12, 13]. The 1/n expansion of critical expo- 
nents |, g, II is based on a similar idea with the only essential difference that 1 jn 
appears as an expansion parameter (at large n) instead of e. Alternatively, it has 
been proposed |ll| to expand the critical exponents in terms of the renormalized 
coupling constant at a fixed dimension d = 3. Later this method has been developed 
by several authors [jl6| , 17, [l^, |l^, ^0], 21]. Apart from the fundamental questions 
concerning the validity of the formal expansion in terms of In k (In k diverges at 
k — ► !) and similar formal expansions which lie in the basis of the theory, a com- 
mon problem for all these methods is that the resulting series for critical exponents 
are divergent (asymptotic), therefore, much efforts have been devoted to develop 



appropriate resummation techniques [18, 19, |2^, 23, 24]. 



In spite of the claims about very accurate values of critical exponents predicted by 
the usual RG theory, we have revealed some serious problems concerning the validity 
of the basic assumptions of this theory. In particular, we have demonstrated that 
the standard RG treatment is contradictory and therefore cannot give correct values 
of critical exponents. Namely, based on a method which is mathematically correct 
and well justified in view of the conventional RG theory, we prove the nonexistence 
of the non-Gaussian fixed point predicted by this theory (Sect. |2|). In Sect. || we 
prove that a correctly treated diagram expansion provides results which essentially 
differ from those of the perturbative (diagrammatic) RG theory. 

Thus, it is worthwhile to search for some alternative analytical methods. One 
of such candidates could be the conformal field theory applied to three-dimensional 
systems [p5|] . Note that in two dimensions this method allows to find the exact crit- 
ical exponents, and even calculate the universal ratios of amplitudes. Some simple, 
but quite plausible models, like the fractal model of critical singularity proposed 
by Tseskis |^ |, also are interesting. We have proposed a novel analytical method 
of determination of critical exponents in the Ginzburg-Landau model (Sec. |||, |5|), 
and have compared the predicted exact values of critical exponents to the results of 
numerical and real experiments (Sec. |(|). 



2 Critical analysis of the perturbative RG method 

Here we consider the Ginzburg-Landau phase transition model within the usual 
renormalization group approach to show that this approach is contradictory. The 
Hamiltonian of this model in the Fourier representation reads 

Y = ^ ( r ° + Ck2 ) ' ^ k ! 2 + uV ~ 1 ^ ^kjVk^ka^-ki-ka-kg , (1) 
k ki,k2,k3 
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where y?k = V^ 1 ^ 2 J exp(— ikx) dx are Fourier components of the scalar order 
parameter field y?(x), T is the temperature, and V is the volume of the system. 
In the RG field theory [|] Hamiltonian ([!]) is renormalized by integration of 
exp(—H/T) over (p^ with A/s < k < A, followed by a certain rescaling procedure 
providing a Hamiltonian corresponding to the initial values of V and A, where A 
is the upper cutoff of the (£> 4 interaction. Due to this procedure, additional terms 
appear in the Hamiltonian (jl|), so that in general the renormalized Hamiltonian 
contains a continuum of parameters. The basic hypothesis of the RG theory in 
d < 4 dimensions is the existence of a non-Gaussian fixed point fi = fx* for the RG 
transformation R s defined in the space of Hamiltonian parameters, i.e., 

Rsfi* = M* . (2) 

The fixed-point values of the Hamiltonian parameters are marked by an asterisk 
(tq, c*, and u*, in particular). Note that fx* is unambiguously defined by fixing the 
values of c* and A. According to the RG theory, the main terms in the renormalized 
Hamiltonian in d = 4 — e dimensions are those contained in (|l|) with Tq and u* of 
the order e, whereas the additional terms are small corrections of order e 2 . 

Consider the Fourier transform G(k, fx) of the two-point correlation (Green's) 
function, corresponding to a point fx. Under the RG transformation R s this function 
transforms as follows []|] 

G(k, n) = s 2 ^ G(sk, R sf i) . (3) 

Let G(k, fj) = G(k, fx) (at k ^ and V — ► oo) be defined within k < A. Since 
Eq. (|3|) holds for any s > 1, we can set s = A/k, which at [i = [i* yields 

G(k, fi*) = a k~ 2+v for k < A , (4) 

where a = A 2 ~' n G(A, (i*) is the amplitude and r] is the universal critical exponent. 
According to the universality hypothesis, the infrared behavior of the Green's func- 
tion is described by the same universal value of rj at any fi on the critical surface 
(with the only requirement that all parameters of Hamiltonian (jl|) are present), i.e., 

G(k, n) = b(ji) k~ 2+r i at k -» , (5) 

where 

b(fx) = ]imk 2 ^G(k,fi). (6) 

According to Eq. (||), which holds for any s = s(k) > 1 and for s = A/k in particular, 
Eq. (|6|) reduces to 

b(n) = lim k 2 ~ r < s{k) 2 - r < G(sk, R s fi) = a, (7) 

if the fixed point fx* = lim R s ll exists. Let us define the function X(k, fx) as 
X(k,fx) = k^G^ik, fx). According to Eqs. (f§), (|), and (|), we have (for k < A) 

X{k,fx*) = - a k^ (8) 
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and 

X(k,n) = -k- 71 + SX(k,fi) , (9) 

where fi belongs to the critical surface, and SX(k,fi) denotes the correction-to- 
scaling term. From (|8|) and (|9|) we obtain the equation 

5X(k, ft* + 6fi) = X(k, fi* + 6fi) - X(k, fi*) , (10) 

where Sfj, = fi — fi* . Since Eq. fllQ| ) is true for any small deviation 5fi satisfying the 
relation 

ff = lim R s {fi* + 8fi) , (11) 

s — >oo 

we choose 5fi such that /x* =^ fi* + 5fi corresponds to the variation of the Hamiltonian 
parameters r^j =>- Tq + 5ro, c* =>■ c* + 5c, and u* u* + e x A, where A is a small 
constant. The values of 5rQ and 5c are choosen to fit the critical surface and to 
meet the condition (|ll]) at fixed c* = 1 and A = 1. In particular, quantity be is 
found 5c = B e 2 + o(e 3 ) with some (small) coefficient B = B(A), to compensate 
the shift in c of the order e 2 due to the renormalization (cf. ||). The formal e- 
expansion of <DT(k, /u), defined by Eq. ([l0[), can be obtained in the usual way from 
the perturbation theory. This yields 

SX(k,fi) = e 2 [Ci(A) + C 2 (A) lnfc] + o(e 3 ) at fc -» , (12) 

where Ci(A) and C*2(A) (C2 7^ 0) are coefficients independent on e. 

It is commonly accepted in the RG field theory to make an expansion like (12), 
obtained from the diagrammatic perturbation theory, to fit an asymptotic expansion 
in k powers, thus determining the critical exponents. In general, such a method is not 
rigorous since, obviously, there exist such functions which do not contribute to the 
asymptotic expansion in k powers at k — ► 0, but give a contribution to the formal e- 
expansion at any fixed k. Besides, the expansion coefficients do not vanish at k — > 0. 
Trivial examples of such functions are e m exp(—ek~ e ) and e m [1 — tanh(e k~ e )] where 
m is integer. Nevertheless, according to the general ideas of the RG theory (not 
based on Eq. (|iol)), in the vicinity of the fixed point the asymptotic expansion 

X(k, fi) = -k~ n + hk^ ^ + b 2 k 2+ °^ + ... (13) 
a 

is valid not only at k — > 0, but within k < A. The latter means that terms of the 
kind e m exp(— ek~ e ) are absent or negligible. Thus, if the fixed point does exist, then 
we can obtain correct e-expansion of <Uf(k, ft) at small k by expanding the term 
bik e+ °^ ) (with b\ = b\ (e, A)) in Eq. (|i~3|) in e powers, and the result must agree 
with ( [12] ) at small A, at least. The latter, however, is impossible since Eq. (12) 
never agree with 

SX(k,fi) =&i(e,A) [l + elnfc + o(e 2 )] (14) 

obtained from ( |T3| ) at k — > 0. Thus, in its very basics the perturbative RG method 
in 4 — e dimensions is contradictory. From this we can conclude that the initial 
assumption about existence of a certain fixed point, predicted by the RG field theory 
in 4 — e dimensions, is not valid. 
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3 A model with quenched randomness 



Here we consider the Ginzburg-Landau phase transition model with 0{n) symmetry 
(i.e., the n-vector model) which includes a quenched randomness, i.e., a random 
temperature disorder. One of the basic ideas of the perturbative RG theory is that 
n may be considered as a continuous parameter and the limit n — ► makes sense 
describing the self-avoiding random walk or statistics of polymers ]|, ||. We have 
proven rigorously that within the diagrammatic perturbation theory the quenched 
randomness does not change the critical exponents at n — > 0, which is in contrast to 
the prediction of the conventional RG theory formulated by means of the Feynman 
diagrams. 

The Hamiltonian of the actually considered model is 



H/T = |[(ro + v^/(x))^ 2 (x) + C (V^(x)) J 



dx 



(15) 



+ uV 1 ^2 ^ i (k 1 )^i(k 2 )'Uk 1 +k 2 <£j(k 3 )<^(-ki -k 2 - k 3 ) 
i j,ki,ka,k3 



which includes a random temperature (or random mass) disorder represented by 
the term y/u /(x) </? 2 (x). For convenience, we call this model the random model. In 
Eq. ( |l5"l) y?(x) is an n-component vector with components <Pi{x) = 
V~ x l 2 J2k<\ ^•i(k)e lkx , depending on the coordinate x, and /(x) = V~ x l 2 / k e lkx 
is a random variable with the Fourier components / k = V~ l l 2 J /(x)e~ lkx <ix. The 
only allowed configurations of the order parameter field </?(x) are those correspond- 
ing to y>i(k) = at k > A. This is the limiting case m — >• oo of the model 
where all configurations are allowed, but Hamiltonian ( |l~5|) is completed by term 
E, k (fc/A) 2m |^(k)| 2 . 

The system is characterized by the two-point correlation function Gj(k) defined 
by the equation 

(^(k)^-(-k)) = 5 id Gi(k) = Si d G(k) . (16) 

It is supposed that the averaging is performed over the y(x) configurations and then 
over the /(x) configurations with a fixed (quenched) Gaussian distribution P({/k}) 
for the set of Fourier components {/ k }, i. e., our random model describes a quenched 
randomness. 

We have proven the following theorem. 

Theorem. In the limit n — > 0, the perturbation expansion of the correlation 



function G(k) in u power series for the random model with the Hamiltonian (15) 
is identical to the perturbation expansion for the corresponding model with the 
Hamiltonian 



H/T 



r <^ 2 (x)+c(V^(x)) 5 



dx 



(17) 



+ uV 1 ^2 ^(ki^iO^^ki+ka yi( k 3)^j(-ki - k 2 - k 3 ) 

i,i,ki,k 2 ,k 3 



where u k = n k - \ (| / k | 2 
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For convenience, we call the model without the term ■ v /u/(x) y? 2 (x) the pure 
model, since this term simulates the effect of random impurities Q. 

Proof of the theorem. According to the rules of the diagram technique, the 
formal expansion for G(k) involves all connected diagrams with two fixed outer 
solid lines. In the case of the pure model, diagrams are constructed of the vertices 
><ww<; ( with factor —uV" 1 ^ related to any zigzag line with wave vector k. 
The solid lines are related to the correlation function in the Gaussian approximation 
G (k) = 1/ (2r + 2ck 2 ). Summation over the components ty?j(k) of the vector y(k) 
yields factor n corresponding to each closed loop of solid lines in the diagrams. 
According to this, the formal perturbation expansion is defined at arbitrary n. In 
the limit n — > 0, all diagrams of G(k) vanish except those which do not contain the 
closed loops. In such a way, for the pure model we obtain the expansion 

G(k) = k -_k + kf^-k + ... . (18) 

In the case of the random model, the diagrams are constructed of the vertices 
>» <C and ••■•< . The factors uV^ 1 (\ /k | 2 ^> correspond to the coupled 

dotted lines and the factors — uV^u^ correspond to the dashed lines. Thus, we 
have 



G(k) = k :k + 



+ ... . (19) 



In the random model, first the correlation function G(k) is calculated at a fixed {/k} 
(which corresponds to connected diagrams where solid lines are coupled, but the 
dotted lines with factors — \fu V^ 1 / 2 f^ are not coupled), performing the averaging 
with the weight f ({/k}) over the configurations of the random variable (i.e., the 
coupling of the dotted lines) afterwards. According to this procedure, the diagrams 
of the random model in general (not only at n — > 0) do not contain parts like , 

> i e tc-> which would appear only if unconnected (i.e., consisting of 

separate parts) diagrams would be considered before the coupling of dotted lines. 
It is evident from Eqs. (18) and (|19|) that all diagrams of the random model are 



obtained from those of the pure model if any of the zigzag lines is replaced either by 
a dashed or by a dotted line, performing summation over all such possibilities. Such 
a method is valid in the limit n — > 0, but not in general. The problem is that, except 
the case n — > 0, the diagrams of the pure model contain parts like , ^CCC j 

, etc. If all the depicted here zigzag lines are replaced by the dotted lines, 
then we obtain diagrams which are not allowed in the random model, as explained 
before. At n — > 0, the only problem is to determine the combinatorial factors for 
the diagrams obtained by the above replacements. For a diagram constructed of M\ 

vertices ^> <C and M2 vertices the combinatorial factor is the number 

of possible different couplings of lines, corresponding to the given topological picture, 
divided by M X \M 2 \. 

Our further consideration is valid also for the diagrams of free energy (at n — > 
represented by the main terms containing single loop of solid lines) and of 2m-point 
correlation function. We define that all diagrams which can be obtained from the 
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i-th diagram (i.e., the diagram of the i-th topology) of the pure model, belong to 
the i-th group. Obviously, all diagrams of the i-th group represent a contribution 
of order u l , where I is the total number of vertices >^^v^<; in the i-th diagram. 
The sum of the diagrams of the i-th group can be found by the following algorithm. 

1. Depict the i-th diagram of pure model in an a priori defined way. 

2. Choose any one replacement of the vertices >*vs<w<; by >=> <C and 

^> <C , and perform the summation over all such possibilities. For any 

specific choice we consider only one of the equivalent M\\M2\ distributions of 

the numbered M\ vertices ^> <C and M2 vertices over the fixed 

numbered positions instead of the summation over all these distributions with 
the weight l/^MxlM^Vj. Thus, at this step the combinatorial factor for any 
specific diagram is determined as the number of possible distributions of lines 
(numbered before coupling) for one fixed location of vertices consistent with 
the picture defined in step 1. 

3. The result of summation in step 2 is divided by the number of independent 
symmetry transformations (including the identical transformation) for the con- 
sidered i-th diagram constructed of vertices ><w^<^ ; since the same (orig- 
inal and transformed) diagrams were counted as different. 

Note that the location of any vertex ^> <C is defined by fixing the position of 

dashed line, the orientation of which is not fixed. According to this, the summation 
over all possible distributions of lines (numbered before coupling) for one fixed loca- 
tion of vertices yields factor 8 Ml 4 M2 / 2 . The i-th diagram of the pure model also can 
be calculated by such an algorithm. In this case we have 8' line distributions, where 
I = Mi + M2/2 is the total number of vertices >*ww<; in the i-th diagram. 
Obviously, the summation of diagrams of the i-th group can be performed with 
factors 8' instead of 8 Ml 4 M2 / 2 , but in this case twice smaller factors must be related 
to the coupled dotted lines. The summation over all possibilities where zigzag lines 
are replaced by dashed lines with factors —uV iik and by dotted lines with factors 
TjuV^ 1 (^\ /k | 2 ^>, obviously, yields a factor uV^ 1 ^— «k + \ (\ /k = —uV" 1 ^ 
corresponding to each zigzag line with wave vector k. Thus, the sum over the dia- 
grams of the i-th group is identical to the i-th diagram of the pure model defined 
by Eq. ([l7]). By this the theorem has proved not only for the two-point correlation 
function, but also for 2m-point correlation function and free energy. 

If, in general, the factor y/u in Eq. ([l^) is replaced by \/vf, where v! is an inde- 
pendent expansion parameter, then our analysis leads to the above relation between 
diagrams for uu^ = uu^ — ^- (\ /k | 2 ^. According to this, at n — > the pure and ran- 
dom models cannot be distinguished within the diagrammatic perturbation theory. 
If, in principle, critical exponents can be determined from the diagram expansions 
at n — > 0, as it is suggested in the usual RG theory, then the same critical exponents 
should be provided for both models at n — > 0. In such a way, we conclude that the 
RG method is not correct because the above condition is violated. As compared to 
our simple treatment of the random model, the RG treatment includes additional 
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Feynman diagrams because the Hamiltonian becomes more complicated after the 
renormalization. However, this does not enable to find the difference between both 
models: the original information, when one starts the perturbative renormalization 
of Hamiltonian (15), is contained in the Feynman diagrams we considered, but the 
renormalization by itself does not create new information about the model. Really, 
by renormalization we merely "forget" some information about the short-wave fluc- 
tuations to make that for the long-wave fluctuations easier accessible. Thus, our 
conclusion remains true. 



4 Equations from reorganized perturbation theory 

As we have already discussed in Sect. ||, it is not a rigorous method to make a 
formal expansion like (^) and to try calculate the critical exponents therefrom. 
We propose another treatment of the diagrammatic perturbation theory for the 
Ginzburg-Landau model defined by Eq. (|i~7|), where uu^ = u^. The basic idea is 
to obtain suitable equations by appropriate grouping of the diagrams. Suitable are 
such equations which allow to find the asymptotic expansions at the critical point 
directly in k power series, but not in terms of the formal parameter In k (as in Eq. 12) 
which diverges at k — * 0. 



4.1 Some fundamental problems of the perturbation theory 

One of the problem which can arise in any perturbation theory is that the pertur- 
bation expansion alone does not define unambiguously the original function. For 
example, all expansion coefficients in the formal expansion of exp(— 1/u) with u 
considered as an expansion parameter at u = +0 are zero, whereas the original 
function is not zero. However, if quantities A and B have the same perturbation 
expansion in a power series of it, then A(u) — B(u) = C(u) holds, where C(u) is 
some function of u with all the expansion coefficients equal to zero. Accordingly, 
the diagram technique principally allows, in the worst case (where C(it) ^ ), to 
find G(k) with an accuracy to some unknown function which tends to zero at u — > 
faster than u l at any I > 0. We have shown in Sect, 
exponents can be obtained neglecting this function. 

Another problem is that the considered formal diagram expansions diverge. Nev- 
ertheless, the following way is possible, which lies in the basis of our diagrammatic 
treatment. First, we build up a perturbation expansion of quantity A which is nec- 
essary to be found. Then we seek such a quantity B (represented by converging 
sums and integrals), the perturbation expansion of which is identical to that we 
have built up. According to the above consideration, we have A(u) = B(u) + C(u), 
where C(it) is an insignificant correction (or zero). In this case it is not necessary 
that the perturbation sum converge if calculated in a straightforward way. Various 
manipulations with diagram blocks appearing in Sec. |^ are defined as constructions 
of corresponding formal (diverging) expressions which, however, provide correct ex- 
pansion in terms of it. Various diagram representations of a given quantity yielding 
the same expansion in a power series of u are defined as equivalent. 



5.1 that true (exact) critical 
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4.2 The diagram notation 

Here we define some diagram notations appearing in Sec. |I[ 

Coupled diagram is denned as any diagram which does not contain uncoupled 
lines, i. e., any line starts from some kink and ends in the same or another kink, where 
"kink" means a merging point of solid and dashed lines of a vertex >» <C . 

The self-energy block — X-@ X~ k denotes the perturbation sum involving all 
connected diagrams of this kind, i. e., the sum of all specific self-energy blocks or 

diagrams which cannot be reduced to a linear chain like ^X-^ — ••■ X"^ 
consisting of two or more blocks. The simplest specific self-energy blocks are 

k T _ k = -2n F- 1 nEGo(q) and kJ...kk = -4V" 1 £ u q G (k - q). 
X i X — q ' q 

Factors corresponding to the lines marked by crosses are omitted. Each case of 
topologically nonequivalent coupling of lines corresponds to one diagram. 

Skeleton diagram is defined as a connected diagram, containing no parts like 
— ^ — , with factors G(k) corresponding to the solid lines. For example, the 
simplest coupled skeleton diagrams are 

O-O =-u V- 1 (£G f (q)j and f^) = -2V" 1 £ E £ Gi(q)Gi(k - q). 

\i,q / j k q 

Note that index % can be removed, replacing J2i with factor n. The simplest skeleton 

diagrams like ^-X-^ ) X k with two outer lines (not coupled diagrams) are the same 
as the above given self-energy diagrams, but with G(k) instead of Go(k). 

Single block is defined as a connected diagram with two outer (broken) dashed 

lines which cannot be reduced to a linear chain like •*• con- 
sisting of two or more blocks. 

The single skeleton block denotes the perturbation sum involving all 

single blocks that belong to skeleton diagrams, i. e., the sum of all specific single 
skeleton blocks. The outer (broken) dashed lines or kinks in this case are marked by 
1 and 2 (in general kinks are not depicted). Factor — V -1 Uk corresponds to the pair 
of these lines. Each case of topologically nonequivalent coupling of lines with respect 
to fixed kinks 1 and 2, considered as nonequivalent, corresponds to one diagram of 

- — . For instance, 

-"d)-~ = 16^ k y- 2 £EE Gi(q)Gi(k - q) n q _ p G i (p)G i (k - p). 

i q p 

Combinatorial factor corresponding to any specific diagram is not given explicitly, 
but is implied in the diagram itself. It can be calculated following the scheme in 
Sec. H (but without the replacements, since now we have only one kind of vertices). 

4.3 Expansion of G(k) in terms of skeleton diagrams 

It is suitably to have a diagram expansion for G(k) where the true correlation 
function G(k) is related to solid lines instead of Go(k). To obtain this, first let us 
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Figure 1: An example of replacement, where the blocks of primary chain inside a diagram 
are replaced by the self-energy blocks, yielding a diagram of set A. The connecting solid 
lines of the primary chain are marked by 1. There is another (not primary) chain, marked 
by 2, which is built into a block of the primary chain. 



consider the quantity S(k) defined by equation 

G (k) = 2^ . (20) 

l-2E(k)G (k) 

It is well known || that terms of the perturbation expansion of S(k) are diagrams of 

the self-energy block — X-^ X" . Note only that the self-energy defined in Ref. || 
corresponds to — 2E(k). 

The desired expansion with G(k) instead of Go(k) is obtained by grouping of 
diagrams involved in the self-energy block. First we consider specific self-energy 

blocks, called the primal diagrams, from which no block of the kind ^ can 

be extracted. Then we consider the set A of diagrams obtained by extending the 
solid lines inside the primal diagrams by adding all possible numbers m £ [0; oo] 
of self-energy blocks to each of these lines. In particular, the diagram of set A, 
depicted on the right hand side of Fig. [I], is obtained from the primal diagram 

ky y-k by adding two (m = 2) self-energy blocks to the inner solid line. 

Obviously, all diagrams of set A are self-energy diagrams, i. e., they cannot be 

split in two blocks like — X-^— @ X k . Besides, any specific self-energy block is 
contained in a perturbation sum represented by one of diagrams of set A. This is 
proved considering primary chains contained inside the specific self-energy blocks. 
A primary chain is defined as a linear chain of specific self-energy blocks which does 
not belong to (i. e., is not built into) a block of some other linear chain inside the 
diagram. If we replace the blocks of primary chains by — ^ — , we obtain a 
diagram of set A, i. e., a diagram consisting of separate linear chains of the new 
blocks, the initial diagram being involved as a particular case (since — ^ — 
involves all specific blocks of such kind). A specific example of such a replacement 
is illustrated in Fig. [l|. Thus, set A contains all specific self-energy blocks. 

Because blocks are distributed independently over separate linear chains, summa- 
tion over all possible lengths A?" of linear chains yields a factor 

J2n Go(k) ^2G (k) x %^-X-^ identical to the perturbation sum of G(k) for 
each of the original solid lines (in a primal diagram) with wave vector k. Thus an 
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equivalent diagram representation of £(k) (Sect. 4.1) is obtained if the perturba- 
tion sum over all diagrams of set A is replaced by the perturbation sum over the 
primal diagrams in which (?o(k) is replaced by G(k). According to the definition in 

Sect. |4.2| , diagrams of the new expansion are skeleton diagrams like . 

Let us assume that there exists a quantity D{G) the perturbation sum of which 
consists of all coupled (containing no outer lines) skeleton diagrams. D{G) have 
to be considered as a function of discrete variables G$(k) (although Gj(k) = G(k) 
holds, the cases with different i formally are nonequivalent) corresponding to the set 

of discrete wave vectors k. All skeleton diagrams of E(k), i. e., those like 
(and these exclusively) can be obtained by breaking a line with wave vector k in a 
coupled diagram of D{G) and removing factor Gr(k) corresponding to this line. This 
procedure is identical to the derivation of the perturbation sum of D(G) with respect 
to Gj(k). Consequently, for any quantity D(G), having the diagram expansion 
represented by coupled skeleton diagrams, the equation 

is true with some function t?(k) providing zero contribution to the expansion in 
power series of u (cf. Sec. |4.1| ). 



4.4 Representation of skeleton diagrams by single skeleton blocks 

In this section the skeleton diagrams are represented by single skeleton blocks. The 

only skeleton diagram among those represented as ^ ^ is O — O ; i- e -> if 

the specific block ^ contains more than one solid line, obtained by coupling of 

two solid lines of the " half-vertex" > , then we can extract the part — ^ — 

by breaking the two solid lines coupling it to the rest of the block. 

Other skeleton diagrams without outer lines, representing the perturbation sum 
of the quantity 

D*(G)=D(G)- O-O (22) 

can be represented in one way, at least, by cyclically coupled specific single skeleton 
blocks. Really, we can break any of the dashed lines in the skeleton diagram of 

D*(G) to obtain a connected diagram like ^ . In this case ^ 

denotes any skeleton diagram with two outer dashed lines. If the obtained diagram 

cannot be split to make two ^ ^ by interrupting one more dashed 

line, we can join the broken dashed line again representing this as a single skeleton 

block ^ / • If the obtained diagram can be split in two such blocks but cannot 
be split in three by interrupting one more dashed line, we can represent this by a 

cycle like ^ ^ with two specific single skeleton blocks, when coupling again, and 
so on. Any such cyclical coupling of specific single skeleton blocks corresponds to 
one of the skeleton diagrams because it is impossible to extract any block of the 
kind — ^ — therefrom. It cannot be extracted from any specific single skeleton 
block according to the definition, as well as in any other way because subblocks 
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of the kind ^ (constructed of odd number of solid lines by coupling) are 

principally impossible (in the opposite case parts like ^ ^ — ••• — ^ 

could be extracted). 



4.5 Summation of the simplest skeleton diagrams 

Consider now a contribution to D(G) (denoted by D^(G)) of dia grams containing 
the simplest specific single skeleton blocks 

£(°>(k)= *-0~ =-2n k y- 1 ^G J (q)G 4 (k-q) (23) 

exclusively. This contribution is not merely a formal sum since at small itk (at 
r > 0) it can be obtained by straightforward summation, i. e., 

r~\ A ,0\ 

d^\g) = a-o + ; o J + ; o ; + + - 

= O-O -^X;in[l-2EW(q)] . (24 ) 
/ q 

Equation (|23) is true in any case if D(°)(G) is defined quantity having this dia- 



gram expansion. The sum (24) is calculated by the following method. To calculate 
the combinatorial factor for a cycle comprised of iV blocks, first we count all possi- 
ble distributions of N numbered vertices over N sequentially numbered fixed sites 
along the cycle with all possible distributions of 4N uncoupled lines, which are then 
coupled in 2N lines. Then, the result is corrected taking into account that couplings 
obtained in such a way contain equivalent ones which differ merely by a diagram 
having been rotated as a whole or (and) transformed by a mirror-symmetry trans- 
formation, as well as in 2^ ways transformed by a mirror-symmetry transformation 

of each — {3 block separately. 2N-2 N such independent transformations exist. 

4.6 Grouping of the skeleton diagrams 

In this section we reorganize the perturbation expansion which is necessary for cal- 
culation of D*(G). 

First, we have got a formal diagram equation performing manipulations similar to 
those in Sect. |4.5| , but including all diagrams of the single skeleton block. A problem 
arises with more complicated diagrams than those considered in Sect. |4.5| because 
of additional symmetry related to different representations of a given diagram by 
specific cyclically coupled single skeleton blocks. In Fig. @, an example of such a 
transformation is shown, which leads to a different representation of the diagram, 
retaining the couplings of lines unchanged. Two representations of a given coupled 
diagram with numbered vertices are considered as different or nonequivalent if they 
cannot be transformed into each other without using a transformation of this kind. 

If we now consider the perturbation sum of cyclically coupled blocks ^ 



with the combinatorial factors calculated as in Sect. 4.5 without account for the 
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Figure 2: An example of transformation leading to a different representation of the diagram. 



above considered additional symmetry, then a specific diagram of the i— th topology 
(or the i— th skeleton diagram) in the perturbation sum will have weight mi, where 
rrii is the number of possible different representations of this diagram by specific 
cyclically coupled single skeleton blocks. This leads to the following formal identity 
of perturbation sums 



1 



1 - 2 x 



(25) 



where Q denotes the i— th coupled skeleton diagram, the apostrophe " / " showing 
that diagram O — O is omitted. 

Number m; can be identified with the number of linear chains consisting of 

dashed lines and blocks of the kind ^ (or a single dashed line) contained 

within diagram Cj, counting all such chains contained inside specific single skeleton 
blocks that comprise the cycle in a given representation and the chain representing 
the cycle itself. It holds because any linear chain contained in a single skeleton block 
can represent a cycle in one of the new representations of the diagram obtained by 
such a transformation as illustrated in Fig. Note that linear chains can be built 
into each other in all possible ways (like in a fractal), and all these chains (of all 
levels) are counted as different. 

We have avoided explicit counting of linear chains as follows. We have proven the 
following: If there exists a quantity £(q, £) which has the perturbation expansion 
^(QjC) = Z)^i( c l)C n S where £,(q) is a specific single skeleton block of the i— th 

i 

topology and n, is the number of linear chains contained in this block, and if there 
is a solution of the equation 

D*(G, = -\ £ In [1 - 2S(q, C)] - C^D*(G, (26) 



with the boundary condition 

D*(G,0) = -±$> [l 



2£(°>( 



(27) 



then quantity D*(G,C), calculated from this equation, has the perturbation expan- 
sion D*(G,() = X)CiC m,_1 - Since Eqs. (^) and (^7j) define the solution unam- 

i 

biguously, it is necessary only to prove the identity of perturbation expansions in 
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the case if D*{D,Q) has the expansion 2~^CjC mi_1 - The identity of diagram ex- 

i 

pansions in (^) corresponding to m; = 1 follows directly from the consideration 
in Sect. [4.5| . The perturbation expansion of the logarithmic term in (^) contains 
diagrams constructed of specific single skeleton blocks Sj(q) supplied by factors ( ni . 



At ( = 1 this expansion is identical to that on both sides of Eq. (25). Now, any 
diagram is supplied with an additional factor which is the product of factors of kind 
Q ni coming from all specific single skeleton blocks of this diagram, yielding the re- 
sulting additional factor Q™ 1 ^ 1 (consistent with the definition of rrij) or, according 
to (|25|), the weight ^ nii ~ 1 of the i-th skeleton diagram. On the other hand, the 

same diagram appears in the perturbation expansion Ci C mi ~ °f D*(G, C) with 

i 

the weight £ m * _i and in the perturbation expansion of (4^D*(G,() with the weight 
(rrii — 1) C mi ~ ■ This, obviously, leads to the identity of perturbation expansions in 



the left-hand side and in the right-hand side of (26), which proves the statement. 
We have defined D*{G) by 

D*(G) = D*{G,1) , (28) 

which, obviously, provides the perturbation expansion of D*(G) defined before. 
Based on Eqs. (|26|) to (^), we can calculate D*(G) without explicit counting of 
the linear chains. 

At the final step of our diagrammatic transformations, we make the summation 
over the linear chains contained in the single skeleton block following the method in 

Sect. |4.3[ , with the only difference that blocks ^ , instead of — ^ — , 

are considered. 

As a result, the new perturbation expansion of X(q, C) is 



S(q, = ? -0- + < - q -<|>" + < 2 { "lii/-- < 29 > 

+ - q -<X> :q - + - q -4> :<! + - q -oxF- }•• 

where the perturbation expansion of the waved line represents the sum over the 
linear chains of all possible lengths N, i. e., the sum of the geometrical progression, 

= f -uq^7[l " 2E(q, C)] = -u^V- 1 ^ [2E(q, C)]^ • (30) 

N>0 

Acording to this procedure, the new diagrams do not contain parts like — ^ 

and/or v/xo^xru , i. e., they are single blocks and skeleton diagrams with respect 
to both solid and waved lines. By this the grouping of diagrams is completed. Such 
a grouping of diagrams is unique in the sense that it allows to analyze all diagrams 
simultaneously, considering asymptotic equations related to critical phenomena. 

It should be noted that the actual grouping of skeleton diagrams cannot be 
extended trivially to include even the simplest y? 6 term. If we include, e. g., the tp 6 
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vertex ^* ^ , subblocks of the kind ^ are possible ( O , in 

particular), therefore, the cyclical coupling of single skeleton blocks ^ (see 

Sect. yield diagrams which are not skeleton diagrams anymore. For example, 

the diagram represents one of the diagrams involved in ^ . 

The cyclical coupling of two such blocks is not a skeleton diagram because it can be 
split in two blocks O O with two outer solid lines. 



4.7 Equations for calculation of the correlation function 

In this section equations are considered from which S(q, C) and G(k) can, in prin- 
ciple, be calculated. Although this quantity is not defined unambiguously by the 
perturbation expansion, we can use any of possible functions having the expan- 
sion Q29|). Only function i?(k) in Eq. ( |2~T| ) can be finally affected by the specific 
choice. We have defined X(q, C) by equations 

u~ p u-p 

S(q,C) = £(°)(q)+ J e'^dh J e^B^, C, *i* 2 ) dt 2 , (31) 



o 



00 ^m+m 



(ml) 

771=1 v ' 

where S^ m ^(q, C) represents the sum of diagrams of the m-th order (m > 0) in 
(p9j), and p is a constant having the value < p < 1/2. Term £(°)(q) in ( |3l|) 
is separated to ensure that the boundary condition S(q, 0) = S^°)(q) is satisfied 
when the first diagram is retained in (p^). Here and in our further considerations 



expansion in the vicinity of the point u = +0 is used. Equations ( |3l| ) and (32) 
yield function S(q, Q) which has the perturbation expansion equivalent (i. e., the 
same, if represented in terms of u) to (|29|), and to S(q, C) = X)5]i(q)C n S since 



u-P 

J e ~H m dt = ml + o (u- pm exp (-u- p )) and B(q,C,*) ~ Cut (for < t < u~ 2p ) 


hold at u — > and < p < 1/2. The latter ensures convergence of integrals ( |3lD 
at u — > 0. This relation holds because the actual expansion parameter in ( |32[ ) at 
u — ► is Cut and this sum converges absolutely, as discussed below. At large m 
the number of terms of the m-th order increases approximately as m 7 6 m m! with 
some constant 7 and these terms behave as — (— b) m where b > 0. According to this 
consideration terms in ( |3~2| ) can be approximated as const ■m 1 (—bCt) m / ml with new 
(6 times larger) value of b. Consequently, the perturbation sum of S(q, C) diverges 
at any values of parameters, whereas ( |32| ) converges absolutely at any given values 
of u, C an d t. Therefore, we can find -B(q, C> t) from (|32] ) with any desired accuracy 
(the larger are C an d t, the larger is the number of terms to be counted). Then we 
can find £(q, C) from (31). 

The Dyson equation for G(k) following from Eqs. (|20|) and (|2l|) is 

1 -ro + ck*-^-+m- (33) 



2G(k) u dGi(k) 
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The same equation with term #(k) neglected has been obtained in Ref. |2j 

Our further analysis is limited to the case u(x) = uS(x) or Uk = u, where u > 0. 
In this case, from p3|), fl2^), and (p8|) we obtain 



l/(2G(k)) = r + c/fc 2 + 2uG + i2(k) + 0(k) , (34) 

where 

G=( 95 2 (x)) = ^ 1 n^G(k), (35) 

k 

fl(k) = "^(kT • (36) 

All terms in Eq. fl34j) are well defined. D*(G, £) and, consequently, -R(k) is defined 
by Eqs. ( |26"[) to (|32|). According to the definition, $(k) is a quantity which have 
to be included to obtain an exact equation which is satisfied by the exact corre- 
lation function G(k) given by the statistical integrals. We know that #(k) does 
not contribute to the formal expansion of G(k) in u power series. It means that 

lim #(k, u)u~ T = holds for any positive r. Our equations have an obvious phys- 
ic— >o 

ical solution S(k,C) S (0) (k), D{G) ~ and G(k) ~ G (k) at r > 

and u — * 0, which agree with the true (exact) G(k). Analytic continuation to ar- 
bitrary vq value is possible if one starts with a finite volume V and consider the 
thermodynamic limit V — * oo afterwards. 

5 Asymptotic solution and critical exponents 
5.1 Determining the critical exponents at u — > 

We will show in this section how the true asymptotic solution at k — * (at the 
critical point T = T c ) can be found from the simplified equations where u tends to 
zero. The right-hand side of ([34|) vanishes at k = and V — > oo at the critical point 
and, thus, we have 

l/(2G(k)) =cfc 2 + J R(k)-i?(0)+7?(k)-i?(0) at T = T c . (37) 

The asymptotic of the correlation function at k — > in this case is G(k) ~ a k~ x 
where a is a constant and A = 2 — 77, r\ being the critical exponent. The case of the 
spatial dimensionality d < 4 is considered. The correlation function is well described 
by G(k) = ak~ x within some critical region k < k cr it(u), where k cr u(u) tends to 
zero at u — > 0, since at u = the Gaussian approximation with A = 2 is the solution 
of (|37|) for any k. Let us define the effective value of r\ at some k = k by 



rj(k,uj = -d (in [fc- 2 Gf -1 (k)]) /d(lnk)\ __ . (38) 
Then, k cr it(u) is defined by 

\r] (k crit (u),u) - rj\ = e , (39) 
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where e is a sufficiently small constant. According to the universality hypothesis, 
we have Urn rj(k,u) = rj for any positive u. We have also lim rj (u r k cr u(u), u) = r/ 

at r > 0, since the critical exponent is determined at k <C k cr it( u ) if it — >• 0, and, 
therefore, corresponds to the asymptotic solution at k —* 0. Based on a non- 
perturbative analysis provided in Appendix B, we have shown that 

lim (u s /k cnt (u)) = , (40) 

holds at large enough s. Quantity k cr it(u) can be related to the region where the cor- 
relation function is well approximated by the asymptotic expansion G(k) =J2^i & Ai 

i 

of any given number of terms. In this case k cr it(u) is defined by the condition that the 

approximation error at k = k cr it{u) corresponds to the variation of the smallest A; by 

some small e. Following the method in Appendix B, the statement ( |40| ) holds in this 

case too, and coefficients ci(u) in the asymptotic expansion 1/G(k) = c i( u ) k 2X ~ Xl 

l 

(where I > and Ao = A) meet the condition lim (u Sl /ci(u)) = at large enough s/. 

According to the discussed property of q, the following is true. If we as- 
sume that, at k = u r k cr i t (u) and u — ► 0, uj(k,u) = #(k) — i?(0) is either com- 
patible with some term q k 2X ~ Xl , or is much larger than all these terms, then 
lim [u T /u (u r k cr it(u) , u)] = holds for large enough values of r. On the other hand, 

the basic property of u(k, u) is that lim (u T / \ u(k, u) I) = oo holds for any r at any 

arbitrarily small, but fixed k (see the end of Sect. [4.7| ). From this and Eq. ( f40| ) we 

conclude the following: if the above assumption is true, then, at u — > 0, | u>(k,u) \ 

decreases faster than any negative power of k if k is increased within some region 

of infinitely small, but larger than u r k cr it(u) values of k. On the other hand, the 

only essential singularity of the correlation function is at k = 0, which means that 

the above assumption leads to unphysical conclusions regarding behavior of G(k), 

calculated from Eq. (|37|), unless the sharp decrease of | u(k,u) \ is compensated by 

the corresponding variation in -R(k) — i?(0). However, our further analysis strongly 

supports the idea that i?(k) — R(0) is a well defined smooth function of k which 

behaves like some power of k at k — > 0. Thus, the compensation is impossible and 

the discussed here assumption is false, i. e., the opposite is true: at k = u r k cr it(u) 

and u — > the term co(k,u) is negligible compared to any of corrections to scaling 

in the asymptotic expansion of 1/G(k). 

In principle, critical exponents can be found by calculating G(k) within k S 

[u r k cr it(u); A] directly from Eq. (]37|) where the term #(k)— 1?(0) is neglected, followed 

by extrapolation of the results to smaller values of k in the form G(k) = J2^i k~ 1 . At 

l 

u — > this yields true (exact) correlation function, since Eq. ( |37| ) is exact (according 
to the definition), term i?(k) — t?(0) is negligible, and the lower marginal value of 
the considered interval is infinitesimal compared to k cr it(u). According to the above 
analysis, this method yields exact critical exponents Aj. 
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5.2 Scaling properties of the main terms at T = T c 



It is impossible to calculate precisely all terms in Eq. (|32|). However, this is an unique 
feature of our reorganized diagram expansion that all terms have common scaling 
properties where the order of diagram does not appear as a relevant parameter. Our 
claims are based on the proof of these scaling properties at d = 2, 3 for separate terms 
and also for the whole sum (|32|) not cutting the series. To simplify the notation, n 
is considered as a fixed parameter not included in the list of arguments. 

We introduce a lower limit of the wave vector k m in to simulate a finite-size effect 
(as if the linear size of the system would be L = 2ir/k m i n ) in calculation of -R(k), 
which ensures the convergence of any k-space integral. Finally, we consider the limit 
kmin — >• in the equation for i?(k). 

Since r\ > holds for d < 4, the term ck 2 in (|37|) gives a small correction which 



is by factor kP smaller than the main term R(k) — R(0). Term "1" in (Pq), (27), and 
(p0|) causes a small correction to R(k) as well. This can be checked easily by a direct 
calculation as in Ref. [ 28 [ , if only the first term in (|29|) is retained. We prove here 



that a selfconsistent solution of our equations can be found where this condition is 
satisfied with account for all terms. 

First, let us consider the dominant behavior of £(q, £) calculated at G(k) = 
ak~ x . A property of the asymptotic solution is such that 

E(q, C) = a 2 x q d - 2X ¥ {A/q, k min /q, (, A, d, u) (41) 

holds, where \& is a function of the given arguments only. Such a solution is possible 



since each term in (f32|) has the scaled form like (|4l]) if ([3C]) is substituted by fl4T|) 
neglecting "1". It becomes obvious after the following manipulations: all the sums 
are replaced by integrals (the standard procedure), V -1 J2k = ^ (27r) _d J d d k, and 



all the wave vectors in (|29| ) or (32) are normalized to the current value of q. This 
produces factor aq~ x corresponding to each of the solid lines, factor a~ 2 q 2X ~ d corre- 
sponding to each of the waved lines, and factor q d corresponding to each integration 
for any diagram of (|29|). As a consequence, the resulting factor is a 2 q d ~ 2X irrespective 
of the order of diagram m (there are 2m + 2 solid lines, m waved lines, and m + 1 in- 
tegrations). Thus, a 2 q d ~ 2X appears as a common prefactor for the whole sum (32), 



and the problem reduces to summation of scaling functions depending merely on 
A/q, k m i n /q, C, A, d, and u, which leads to (fi"l|). In this case arguments A/q and 
kmin/q. represent the upper and the lower limits of integration for normalized wave 
vectors related to the solid lines. 

Formally we could allow other kind of solutions, but the true result must coincide 
with the asymptotic expansion E(q, C) = E m C m £( m )(q) at C -> (consistent with 
the boundary condition S(q, 0) = S"(q)), which shows that ( |4l| ) is the only possi- 
bility. If the correction term " 1" in Eq. (|30|) is omitted, then any of the expansion 
coefficients has the scaled form with factor a 2 q d ~ 2X multiplied by some function of 
A/q, k m i n lq, A, d, and u. This is proved by induction over m: it holds at m = 0; if 
it holds for terms up to the m-th order, then it holds for terms of the (m + l)-th 
order, calculated from diagrams with no more than m + 1 waved lines (expanded in 
terms of Q). Thus, the asymptotic solution at Q — > is unambiguous and has the 
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scaled form (|4l|), which is a property of the solution as well at finite £. 

It is purposeful to seek a solution for dE(q, Q/dGiQs) in form 
V- 1 ak- x Y (q/fc, k/fc, A/fc, fc min /fc, C, A, d, u) using Eqs. @ and (||), where S(q, C) 
is considered as a quantity which is already known. The latter can be represented 
in the form a 2 k d ~ 2X \& (q/k, A/k, k m i n /k,(, X,d,u) obtained from ( |4l| ) by changing 
variables q — > k ■ (q/k). This yields a selfconsistent equation for the unknown 
function Y. The arguments of this function only are contained therein, since all 
terms <9£( m )(q, C)/<9G;(k) in the sum for dB(q,(,t)/dGi(k), obtained from (pi), 



have the form V 1 ak x multiplied by some function of these arguments. In a 
similar way, we can find from the closed equations (|26j) and a selfconsistent 
solution dD*(G,()/dGi(k) = a~ l k x 4>(A/k,k min /k,C,\,d,u). The analysis of the 
asymptotic solution at £ — ► shows that discussed above are the true solutions 
meeting fl27[) . 

Consider now the limit k m i n — ► 0. The nonzero lower limit of integration has 
been introduced to avoid the divergence of k-space integrals in (^) considering the 



formal expansion in terms of £. If the selfconsistent solution of Eqs. ( |3l[) and (32) is 
considered, then convergence of integrals at k m i n = is ensured since these equations 
provide the solution of physical problem where the existence of thermodynamic limit 
{kmin — ► 0) is doubtless, and this conclusion agree with the formal analysis of our 
equations at fc m j n = 0: S(q, () diverges at q — » 0, but this term appears in the 
denominator of ([30]) , which ensures the convergence of k-space integrals in ( f2~9| ) . A 
selfconsistent solution with exponentially diverging E(q, £) at q — > is not possible, 
which ensures the convergence in (|26|). Thus, we can set k m i n = 0, which at £ = 1 
leads directly to the relation ak~ x R(k) = <p* (A/k, A, d, u), where (f>* is a function 
exclusively of the given arguments. Since R(0) is a constant, the function eft* can be 
represented as <f>* (A/k,X,d, u) = aA- x R(0) ■ (A/k) x + 1/ [20 (A/k,X,d, u)]. We get 



G(k) = ak- x cP(A/k,X,d,u) (42) 



by substituting this into (I 



According to the consideration in Sect. 5T| , the true value of rj can be found from 



our equations in the limit u — > 0. In this case behavior of G(k) within some region 
k ~ u r k cr it(u) and extrapolation to smaller values of k have to be considered. It 
follows from the definition of k cr n(u) that a stable solution in the form G(k) ~ ak~ x 
does exist at k ~ u r k cr it(u). It means that cf) does not depend on A/k (or k) within 
this region at u — > 0. Thus, the value of argument A/k in ( f42| ) can be replaced 
by Au~ r k~^ it (u). On the other hand, universal positive value of rj is obtained from 



Eq. ( |37|) at u — > 0, which means that a finite limit lirn^^Au r k crit (u), X,d,u 
b(X, d) exists. So, we have an asymptotic equation 

G(k) = a k~ x = a k~ x x b(X, d) (43) 

from which, in principle, the universal value of A can be found ( i. e., b(X,d) = 1). 
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5.3 Scaling properties at T = T c including correction terms 

The analysis of the previous section can be extended by including corrections to 
scaling. As it was mentioned, correction e(k) ~ kP is introduced by ck 2 in ([3?]). As 
regards term "1" in (^), (p7[), and (p0|), it yields (performing the same analysis as for 
the main terms) a correction 5(k) ~ k 2 ^ d in (37) corresponding to the contribution 
of the correction term (with the factor q 2X ~ d ) in equation 

1 - 2S(q, C) ^ -2S(q, Q X (l - q 2X - d / ha 2 * (A/q, k min /q, £ A, d, u)] ) . (44) 



Here, and in equations (|29|) to (|32j) q, A/q, and k m i n /q are considered as independent 
variables, considering the limit q —* at any given A/g, and k m i n /q. Similarly, in 
equations for R(k) and G(k) /c, A/A;, and k m i n /k are considered as independent 
variables. Finally, we set k m i n /k — > to fit the thermodynamic limit, and then we 
consider the limit u — > with the simultaneous tending of A/k to infinity (as in 
Sect. |5.2j ) to get the asymptotic solution at a fixed A. The solution can be expanded 
in terms of e(k) and 5(k) to yield an asymptotic expansion at k —* 

G(k) = J>fc- A < (45) 



where 



\l = \ — ni ■ r) — mi ■ (2A — d) (46) 



with n/, = 0, 1, 2, ... corresponding to the correction of order e n '5 mi , 6/ being 
expansion coefficients. 

If the right hand side of Eq. (|36[) is substituted by (|45|) , we get 



fl(k) = X! 0* ( A A, A, d, u) A: 2A - A ' (47) 
i 

where (f>* (A/k, A, d, u) are functions of the given arguments. They depend also on 
the set of coefficients bi . This is obtained by finding selfconsistent solutions similarly 
as in Sect. |5.2| , with the only difference that any given number of corrections to the 
scaling is included. 

Since R(0) is constant, <j>* behaves as const ■ (A/k) 2X ~ Xl at k — > 0, and herefrom 

it follows that R(0) can be represented as R(0) = J^Ri A 2A_Ai with Ri depen- 

l 

denging on A, d, u, and coefficients b[. Based on the same logic as in Sect. 15.21, 
functions <j>i (A/k, A, d, u) = ^ (A/k, A, d, u) — Ri (A/k) 1 may be replaced by 
4>i (Au~ r k~ r \ t (u) , \,d,uj at u — > to obtain the asymptotic solution (45) at a fixed 
A. This replacement is justified by the following argument. If the amplitude a (= bo) 
is considered as known (fixed) quantity and k ~ u r k cr it(u) — > 0, then <pi (A/k, A, d, u) 
tends to some function 4>i(X, d, u) and 4>i (A/k, A, d, u)—(pi(X, d, u) tends to zero faster 
than any positive power of k, since only in this case exponents in the right-hand 
side of ( |37| ) are not affected by short-wave fluctuations and are the same as those 
in the left-hand side of the equation. 

To show this precisely, we prove the following statement: if in the considered 
limit (pi (A/k, A, d, u) tends to 4n(X,d,u), then the tending is faster than (A/k)~° , 
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where a is any finite and positive constant. Really, if G(k) in the right hand side of 
Eq. © is replaced by G(k) - 5G(k), where <5G(k) = G(k) • {k/A') m / [1 + (k/A') m ], 
A' < A, and m — ► oo, then this is equivalent to the shift of the upper limit of wave 
vector magnitude from A to A' at a fixed amplitude a. On the other hand, the 
quantity 5G(k) can be treated as any other correction term. If for arbitrary (i-th) 
correction cpi tends to fy, then at k ~ u r k cr i t (u) — > the quantity <5G(k) produces 
in the right hand side of Eq. (|37| ) a correction term of order k x+m where m — > oo. 
Thus, the shift of A produces a correction smaller than k a at any finite a, which 
proves the statement. 

According to physical arguments, the condition that (j>i tends to (fii, obviously, is 
satisfied. It means that the main contribution to <pi (A/A;, A, d, u) in the equation for 
1/G(k) at k —* is provided by the integration over small wave vectors (therefore the 
result is almost independent on A), i. e., the critical behavior is governed by the long- 
wave fluctuations. Besides, this condition means that the asymptotic solution in 
approximation G(k) = J2h=o b k~ Xl , including m correction-to-scaling terms (m=0, 
1, 2, etc.), is stable with respect to (i. e., is not changed by) higher order corrections, 
which is reasonable and expected. 

As regards corrections, we allow a possibility that any of corresponding expansion 
coefficients can be zero, our analysis is correct in this case. 



5.4 Asymptotic solution at T — > T c 

In this section we have extended our scaling analysis to describe the critical behavior 
of the model when approaching the critical point from higher temperatures, i. e., at 
positive A = T — T c — ► 0. 

The critical exponents cannot be affected by short-wave fluctuations. According 
to this idea, the contribution of sufficiently large k may be neglected in equations 
for 1/G(0) and [1/G(k)] - [1/G(0)] 

= ^' A +yE t G ( k ) " G *( k )] + R (°) ~ R *(°) » ( 48 ) 
' R(k) - R(0) + ck 2 (49) 



2G(k) 2G(0) 

obtained from Eq. ( |34[ ) omitting the irrelevant correction #(k) and assuming that 
tq(T) = ro(T c ) + (drg/dT) • A is the only parameter in ([17]) which depends on 
temperature. In Eq. fl4S|), G*(k) is the value of G(k) at T = T c and R*(0) is the 
value of i?(0) calculated at G(k) = G*(k). This equation represents the condition 
that 1/G(0) vanishes at T = T c . 

Considering the solution at T = T c , we have concluded that terms in the right- 
hand side of equation for l/(2G(k)), calculated at a fixed G(k) (i. e., at fixed am- 
plitudes bi in (gy)), are not sensitive to a variation in the upper limit of the wave 
vector magnitude A at u — > if A/k ~ u~ r k~ r \ t (u) holds, where r is any positive 
constant. In this case unambiguous solution insensitive to the short-wave fluctua- 
tions is obtained, based on the asymptotic expansion (|45|), if the main amplitude 
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bo is considered as a known (fixed) quantity. The asymptotic solutions at T = T c 
and T > T c join at k ~ l/£, where £ ~ A _I/ is the correlation length, therefore, the 
solution of Eqs. ( |48| ) and ( |49| ) at T — > T c for fc ~ l/£ is insensitive to variation of A 
within some region A ~ £ -1, u _r /i;~^(u) ( a * u ~^ 0), if calculations are performed at 
fixed G*(k) (i. e., at fixed bo). It is supposed that k cr i t {u) is determined at a fixed 
upper integration limit A', and A is smaller than, but comparable with A'. The latter 
condition is satisfied for values of A which are smaller than, but comparable with 
B u r / v A cr it(u), where B is appropriate constant and A cr i t (u) ~ [A(u)k cr u (u)] 1 ^ is 
the width of the critical region inside of which the correlation function is described 
by G(k) = £ A g(k£) with the relative error not exceeding some small given value. 
According to our definition £ ~ A(u)A~ v holds, and the above relation for A cr n[u) 
is true since the width of the critical region for £ _1 is proportional to k cr n(u) due 
to the joining of asymptotic solutions at k ~ l/£. We can conclude from this 
discussion that the region k > C(u)/£ with C(u) = u~ r k~ r ] t (u) corresponds to neg- 
ligible short-wave fluctuations if the solution inside the asymptotic region k ~ l/£ 
is considered at u —* and A ~ B u r ^ u A cr it(u). In such a way, we may neglect 
the short-wave fluctuations by formally setting A = C{u)A u (i. e., G(k) = and 
G*(k) = at k> C(u)A u ) in the right-hand side of Eqs. © and © where G*(k) 
has been calculated before this procedure at the true (constant) upper limit A', and 
C(u) = C{u)/A{u). According to the above discussion, this method provides correct 
correlation function at k ~ l/£, which means that it yields true critical exponents. 
In this case the asymptotic solution can be found in the form 

G(k) = ]TA-^(kA-<0, (50) 

l>0 



where exponents 7; are related to those given by fllq ) in the way predicted by the 
scaling hypothesis, i. e., gi(y) behave like gi(y) ^ hy~ Xl at y — > oo to yield (|45| ) at 
A — > 0. Thus, the exponents in (|50|) are 



7 = \v = {2-r 1 )v, (51) 

7; = ni <5i + mi 5 2 , (52) 

where 

5 1 = 2zy-7>0, (53) 

5 2 = 2 7 - dv > . (54) 

Here the main term is given by I = and 7 is the susceptibility exponent. The l-th 
term with I > represents the correction of order e ni 5 m ', where e(A) = A^ 1 and 
5(A) = A fo . 

This result is obtained by a scaling analysis similar to that we have made at 
T = T c . The only difference is that wave vectors are normalized to A", but not 
to the current value of q (in equation for S(q, ()) or k (in equation for G(k)). 
Besides, the lower limit of integration may be set k m i n = from the very beginning 
since there is no singularity at k = 0. In such a way, retaining only the main 
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term, we prove the following scaled form for the relevant quantities calculated at 
A = C(u)A v : S(k,C) = A- 2 T+ d ^(k',C), 9E(q,C)/5Gi(k) = V^Ytf ,U , Q, 
dD*(G,()/dGi(k) = A^r(k',C), and R(k) = ATr(k'), where k' = kA~ u , and at 
any fixed u, d, n, and <7o(kA _y ) the scaling functions depend on the given arguments 
only. Including corrections produced by ck 2 in Eq. fl34|) (correction of order A* 51 ) 
and "1" in Eq. fl3C| ) (correction of order A* 52 ), we obtain the asymptotic expansion 
in the scaled form 

R(k) = Y / A^r l (k') (55) 

l>0 

with exponents 7; defined by Eq. (52). If Eqs. ( f48[ ) and (|^) are substituted by (55) 
at A = G(u)A u , then we obtain a solution where no other correction exponents 
appear. 

5.5 Possible values of critical exponents 

By using the equations and scaling relations obtained in previous sections, here we 
derive our central result - the set of possible values for critical exponents. 

It is reasonable to assume that 7 > 1, which leads to the conclusion that 
(dro/dT) ■ A in Eq. fl4S|), where A = C(u)A u , is compensated by one of the terms 
coming from the asymptotic expansion 

2u V- 1 G ( k ) ~ G *( k ) = Z) B l A 7 ~ 52+7i . (56) 

It means that 

7m = 1 - 7 + 5 2 (57) 

holds at some m > 0. Condition j m < 82 follows herefrom, since 7 > 1. Thus, with 
account for (|^) we have j m = rnS\, and Eq. ([57]) becomes 

mSi-5 2 + j = l. (58) 

We need one more relation to determine the values of critical exponents. We obtain 
this relation from calculation of specific heat Cy assuming the well known hyper- 
scaling hypothesis 

a + dv = 2 . (59) 

The singular part of the specific heat behaves like A~ a and it can be related to the 
singular part of G, i. e., to G — G* (where G* is the value of G at T = T c ), as follows 



r 9 
v(X dA 



V- 1 £ (G(k)-G'(k)) 

fc<A' 



(60) 



where A' is constant. The latter relation follows from thermodynamics, taking into 
account that 

d fF\ dlnZ „ / 2 



dro \T J dr{ 







V(<p 2 (x)) = VG (61) 
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holds where F is the free energy and Z is the statistical sum, consistent with the defi- 
nition of Hamiltonian ( |i~7|) where tq is the only parameter depending on temperature 
and the dependence is linear. 

According to the universality hypothesis, critical exponents do not depend on 
the coupling constant u, therefore the exponents obtained by our method at u —* 
can be used in calculation of the singular part of Cy at a finite u. Coefficients B\ 



in (56) can be changed, not changing G(k) at k ~ l/£, if the solution at A = A' = 
const instead of the formal solution at A = C(u)A u is considered. For instance, a 
variation in Bi due to the change in G(k) at k ~ C(u)A u 3> 1/^ (at u — > 0) can 
be compensated by a contribution coming from R(0) — R* (0) due to the integration 
over G(u)A u < k < A'. Thus, the true values of the amplitudes in the expansion of 
G — G* are unknown, and we allow all the possibilities. 

Thus, consider a finite u. Our basic idea is that the contribution to (|60| ) pro- 
vided by the summation over k > C/£ or k > CA V at C — > oo (we consider the 
limit A — > at a given C, which then is tended to infinity) with the true corre- 
lation functions G(k) and G*(k) (calculated at A = A') cannot change the critical 
exponent a, because the opposite would mean a violation of scaling relations for 
critical exponents. However, we allow that a logarithmic correction can be caused 
by this contribution. Thus, we can replace the summation limit k < A' in (|6^) by 
k < CA V . The first non-vanishing (i. e., having nonzero amplitude at C = oo) sin- 
gular term in the resulting asymptotic expansion represents the leading singularity 
of Cy. Formally, a constant contribution also is defined as singular (with a = 0) in 
the case if there is a jump of Cy at T = T c from one constant value to another, or 
if a refined analysis reveals logarithmic singularity. Thus, according to (|45|), (|50|), 
(jjTp, and (|60[), all possible values of a are given by 

a = 7 m - 7i (62) 

where i > is integer. In this case gi (y)-bi y~ x * (cf. Eqs. © and @) tends to zero 
at y — > oo and the tending is faster than y~ a where a < d, since the opposite would 
mean that the critical exponent a is changed due to the contribution of k > CA U . 
In the marginal case when gi(y) — b{y~ Xi ~ y~ d holds this contribution yields a 
logarithmic correction, i. e., Cy ~ A~ a In A. In such a way our theory provides an 
explanation of the known logarithmic singularity of the specific heat at n = 1 and 
d = 2. 

We have restricted our analysis to 7 > 1 and a > 1 — 7, which is very reasonable 
assumption in view of the known results. From ( |57|) we obtain 1 — 7 = 7 m — 62, 



which, in this case, yields a = j m — ji > j m — £2 or 7^ < 62- Then, combining (58) 
with fl62| ) and (|59|), with account for definitions (|53|), (|54|), and (|52|), the following 
set of possible values for critical exponents is obtained 

= d+2j + 4m _ = 2(1 + m)+j 
7 d(l + m + j) - 2j ' " d(l + m + j) - 2j ' 1 ' 

where m may have a natural value starting with 1 and j is integer equal or larger 
than —771. This result is obtained by proving all the relevant scaling properties not 
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cutting the perturbation serie s ( | 32|) , based on equations which allow to find the 
exact critical exponents (Sec. |5.1[) . Thus, according to these arguments, Eq. (63) 
represents all possible values for the exact critical exponents at d = 2 and d = 3 
in the cases where the second-order phase transition with spontaneous long-range 
ordering takes place. 

Our analysis is not valid at d > 4 since 8\ and 62 (Eqs. ( |53"| ) and (|54|)) are positive, 
i. e., A Sl and A 52 are small corrections, merely at d < 4. Besides, the analysis 
in Appendix B is not true at d = 4. Solutions at natural n (the dimensionality 
of the order parameter) only have a meaning, since our method is strongly based 
on the proof that in the relevant asymptotical region at u — > the solution of 
our equations agree with the exact correlation function defined by the statistical 
integrals. Although the formal perturbation expansion exists at arbitrary n, such a 
method of proof would be meaningless at a not natural n, since the exact correlation 
function is not defined in this case, and we cannot guarantee that a formal solution 
at arbitrary n has all the correct properties (e. g., existence of the secon-order phase 
transition and scaling relation (|59|)) which have been assumed to derive Eq. (|63|). 
Our predictions do not refer also to the case n = 0. This case is exceptional in 
view of our analysis, since the term 2X(q, £) in the denominator in Eq. (|30|) which 
appears as the main term at n > 1 vanish at n = 0. 

In general, different values of j and m can correspond to different (natural) n, 
i. e., j = j(n) and m = m{n). It is easy to verify that at j = and m = 3 
Eq. ([||) reproduces the known exact results in two dimensions. The known exact 
exponents for the spherical model jjj [n = 00) are obtained at j(n)/m(n) — ► 00. 
Although the derivations are true for d < 4, Eq. (|63] ) provides correct result v = 1/2 
and 7 = 1 also at d = 4. It is reasonable to consider d as a continuous parameter. 
This leads to the conclusion that m = 3 and j = are the correct values for the 
case n = 1 not only at d = 2, but also at d = 3. In the latter case we have 7 = 5/4 
and v = 2/3. The nearest values of 7 and ^ provided by Eq. (|63|), e. g., at j = 1 
and m = 3 or at j = 1 and m = 4 are then the most probable candidates for the 
case n = 2. It is interesting to note that our prediction for the singularity of specific 
heat a = for the Ising model (n = 1) agree with that made by Tseskis [^6| based 
on a fractal model. 



6 Comparison with Monte— Carlo and experimental 
results and discussion 

It is commonly believed that all more or less correct Monte Carlo (MC) simulations 
confirm the values of critical exponents obtained from the perturbation expansions 
based on the renormalization group. This is not true. We have found that some 
kind of MC simulations at the critical point, namely, the MC simulations of fractal 
configurations of Ising model |2^] and the MC simulations of the energy density |3(| 
for the XY model in reality do not confirm the results of the RG theory, but provide 
the values of critical exponents which are very close to those we predicted. 

The MC simulations of Ref. [^] allows to determine the fractal dimensionality 
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Figure 3: Fractal dimensionality D of the three dimensional Ising model at the critical 
point simulated by Monte Carlo method (MCS means Monte Carlo steps). The upper and 
lower dashed lines indicate the theoretical values expected from the known and from our 
critical exponents, respectively. 



D (the largest cluster in the relevant configuration has the volume L D where L 
denotes the linear size of the system) which is related to the critical exponents by 
7 = v{2D — d) or, which is the same, r] = 2 — ^/u = d + 2 — 2D. In our opinion, this 
method is better than other more convenient simulation methods, since it provides 
the value of rj as a result of direct simulation, i. e., there are no fitting parameters. 
Besides, the result is relatively insensitive to the precise value of the critical coupling 
(temperature). In Fig. || we have shown the average values of D (the averaging is is 
made over the MC steps from 1 to 10 (except the initial point), from 11 to 20, and 
so on) calculated from the MC data of Ref. (2t| by measuring deviation from the line 
D = 2.48 in Fig. 8 (of Ref. p9]). If properly treated, these simulation data confirm 
the value of r] about 1/8 (or D = 2.4375) consistent with our prediction 7 = 5/4 
and v = 2/3, as it is evident from Fig. ||. The value D = 2.46 ± 0.01 reported in 
Ref. [^] seems to be determined from the upper MC points (Fig. 8 in Ref. [^] ) only 
which are closer to the known theoretical prediction D = 2.48. 



As regards the MC simulations of the energy density E of XY model [ 30 1 at 
the critical point, the true picture can be reconstructed from the simulated values 
listed in Tab. I of Ref. |pp| . Since all the values of E are of comparable accuracy, it 
is purposeful to use the least-squares method to find the optimum value of \jv by 
fitting the MC data to the prediction of the finite-size scaling theory 



E(L) = E + E\L 



--d 



(64) 



where E{L) is the energy density at the critical temperature T\ depending on the 
linear size of the system L. The standard deviation of the simulated data points 
from the analytical curve (|64]) can be easily calculated for any given value of \jv 
with the parameters Eq and E\ corresponding to the least-squares fit. The result 
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Figure 4: The standard deviation a vs the value of 1/v used in the least-squares fit of the 
finite-size scaling curve to the simulated results including 11 data points (solid curve) and 
9 data points (dashed curve). Minimum of the solid curve, shown by a vertical dashed line, 
corresponds to the best fit 1/v = 1.4457 which is close to our theoretical value 13/9 indicated 
by a vertical dotted line. Other vertical dashed line indicates the value 1.487 proposed by 
authors of Ref . |3C(| . 



is shown in Fig. ||. The thick solid curve is calculated including all 11 data points 
(L=W, 15, 20, 25, 30, 35, 40, 45, 50, 60, 80), whereas the dashed line - including 
9 data points (except L=10, 15) used for the fitting in Ref. |j(J. Minimum of the 
solid curve, shown by a vertical dashed line, corresponds to the best fit 1/v = 1.4457 
which comes very close to our theoretical value 13/9 ( provided by ( |63| ) at j = 1 
and m = 3) indicated by a vertical dotted line. We have estimated the statistical 
error of this MC result about ±0.007 by comparing the best fits for several random 
data sets. Different data sets have been generated from the original one by omitting 
some data points with 10 < L < 80. We have found it unreasonable to omit the 



data points with two smaller sizes, as it has been proposed in Ref. [30], since the 
result in this case becomes very poorly defined, i. e., the dashed curve in Fig. |] 
has a very broad minimum. Besides, there is no reason to omit the smallest sizes, 
since the analytical curve (^4|) excellently fit all the data points and the standard 
deviation for 11 data points is even smaller than that for 9 data points (see Fig. Q). 
The possible systematical error due to the inaccuracy in the critical temperature 
T\ = 2.2017±0.0005 (the error bars are taken from the source of this estimation |31||) 
used in the simulations |3(| has been evaluated ±0.017 by comparing the simulation 
results at T\ values 2.2012, 2.2017, and 2.2022. In this case the values of the energy 
density at a slightly shifted temperature have been calculated from the specific heat 
data given in Tab. I of Ref. |30|. In such a way, our final estimate from the original 
MC data of Ref. [^] is 1/v = 1.446 ±0.025 in a good agreement with our theoretical 
value 13/9 = 1.444... and in a clear disagreement with the usual (RG) prediction 
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Figure 5: Our fits to the original MC data of Ref. for the energy density (left) and for 
the specific heat (right) depending on the linear size of the system L. 



about 1.492. One can only wonder where the value 1.487 proposed in Ref. [ 30 1 comes 
from. It does not correspond neither to the best fit for 11 data points nor to that 
for 9 data points, as it is evident from Fig. |j. The values of 1/u and ajv cannot be 
determined independently from the discussed here energy density data. One of them 



have to be calculated from the scaling relation ajv + d = 2/ v. If authors of Ref. [3C] 
were able to determine \ jv with ±0.081 accuracy, then they should be able to find 
ajv with ±0.162 accuracy. In this aspect, the estimate ajv = —0.0258 ± 0.0075 
given by the authors looks more than strange. 

In Fig. | we have shown our fits to the MC data for the energy density E(L) = 
2.0108 - 2.0286 L" 14 / 9 and for the specific heat c{L) = 7.360 - 6.990 L~ x l 9 . They 
do not look worse than those in Ref. |3(|, but our fit for c(L) seems to be better. 

One believes that the value of critical exponent v about 0.67, predicted by the RG 
theory at n = 2, is well confirmed by very accurate measurements of the superfluid 
fraction p s /p = y in 4 He. This is not true, since in reality these experiments |3^] 
provide a good evidence that the effective critical exponent u e ff(t) = <9(lny)/<9(lni) 
remarkably increases when the reduced temperature t = (T\ — T)/T\ (where T\ is 
the critical temperature) is decreased below 10 -5 . According to Ref. p2| , p s /p is 
given by 

p s /p = y(t) = fco(l + fat){l + D p t A )t< x (1 + 5(t)) , (65) 

where ko, k±, D p , and £ are the fitting parameters, A = 0.5 is supposed to be the 
correction-to scaling exponent, and S(t) is the measured relative deviation from the 
expected theoretical expression obtained by setting 5(t) = 0. The percent deviation 



discussed in Ref. [32] is 100 times 5(t). From Eq. (|65|) we obtain 



, ht AD p t A 1 dS(t) . , 

For the values of t as small as t < 10~ 5 and for 5(t) <C 1 Eq. (|6"fj|) with the fitting 



parameters ( = 0.6705, k = 2.38, fa = -1.74, and D p = 0.396 used in Ref. || 
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Figure 6: Percent deviation of the experimental p s /p data |}2| from the expected theoretical 
relation (Eq. (J6q) at 8 = 0). The stright line shows the slope of this plot at the value of t, 
equal to 5 • 10 , indicated by a vertical dashed line. 



reduces to 



Veff(t)~C + d6(t)/dQnt) 



(67) 



The second term in this equation is proportional to the slope of the percent deviation 
plot 100 6(t) vs hit or lgi (the decimal logarithm) in Figs. 2 and 3 of Ref. Q. We 
have read the experimental data from Fig. 2 in Ref. within the region t < 10~ 4 
and have depicted them in Fig. ^. Almost all the data points with a reasonable 
accuracy fit the smooth curve S(t) vs Igt (dashed line) having a maximum at about 
lgt = —5.5. It means that d 2 5(t)/d(lnt) 2 is negative within some region around the 
maximum, i. e., according to (|67|) the effective critical exponent v e ff(t) increases if 
t is decreased. We have roughly estimated and have shown by stright line the slope 
of this curve at t = t* = 5 • 10~ 7 (t* value is indicated in Fig. || by vertical dashed 
line). From this we obtain d5(t)/d(lnt) ~ 0.025. This result depends on the shift 
in the experimentally determined T\ value. To obtain a more reliable estimate, we 
have performed the same manipulations with the data depicted in Fig. 3 of Ref. |32| | 
corresponding to T\ shifted by ±2§nK, and have obtained the values of dS(t) / <9(ln t) 
about 0.03 and 0.015, respectively. Our final result 0.0233±0.0083 for this derivative 
at t = t* has been obtained by averaging over the three above discussed estimates 
(0.015, 0.025, and 0.03) with the error bars large enough to include all these values. 
According to this, from Eq. @ with ( = 0.6705 we obtain u eff (t*) = 0.694 ±0.009 
which, again, is in a good agreement with the value v = 9/13 ~ 0.6923 provided by 
Eq. (|63|) at j = 1 and m = 3 and in a disagreement with the RG predictions. 
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7 Conclusions 



We have proposed a novel method in critical phenomena (Sect. |], |5|) which is based 
on the grouping of Feynman diagrams in ip 4 model with O(re) symmetry. As a 
result, equations for calculation of the two-point correlation function have been 
obtained containing an infinite, but converging perturbation sum. It has been shown 
that these equations allow, in principle, to find the exact critical exponents. In 
distinction to the usual renormalization group approach, our predictions are based 
not on evaluation of some of the first terms in the perturbation expansion, but on 
the proof of relevant scaling properties for the whole sum, which is possible due to 
the actually proposed reorganization of the perturbation theory. 

Based on this scaling analysis, we have derived a set of possible values for the 
exact critical exponents (|63|). A disagreement with the actually accepted values of 
the critical exponents in three dimensions has been revealed. However, we argue 
that our result is correct since, in distinction to the usual treatment critised in 
Sec. [2| and ||, our method is faultless from the mathematical point of view. Some 
assumptions have been made, but they look innocent and have been well motivated. 
Besides, our method, being equally valid in two and three dimensions, reproduces 
the known exact critical exponents at d = 2. 

A comparison of results has been made in Sect. ^, showing that in some cases, at 
least, our predictions are in accurate agreement with properly treated MC simulation 
data as well as with experiments. More comparison with MC data is in progress. 

In summary, we conclude the following. 

1. The conventional method of the perturbative RG theory is contradictory and, 
therefore, cannot give correct values of critical exponents. 

2. Our equations, derived in Sec. || by grouping of the Feynman diagrams, allow to 
find the exact critical exponents for ip 4 model with 0(n) symmetry (re > 1) by 
proving all the relevant scaling properties of the asymptotic solutions at T = T c 
and T — > T c , not cutting the perturbation series. These scaling properties have 
been proven in Sec. 0. 

3. In the cases of the second-order phase transition with spontaneous long-range 
ordering, all possible values for the exact critical exponents at d = 2, 3 and 
n = 1,2,3, etc. are given by Eq. (|63|). 

4. At m = 3 and j = our result (|63|) reproduces the known exact critical 
exponents in two dimensions at n = 1 (7 = 7/4 and v = 1). Based on the 
idea that d may be considered as a continuous parameter, we conclude that 
7 = 5/4 and u = 2/3 are the true (exact) values at d = 3 and n = 1. 

5. The comparison with Monte-Carlo data in Sec. || well confirms the hypothesis 
that 7 = 17/13 and v = 9/13 (corresponding to m = 3, j = 1) are the true 
values of the critical exponents at d = 3 and n = 2. 
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Appendix A. Large— order estimation 

Let us consider coupled diagrams, constructed of vertices ><'vv\*< , which are 
skeleton diagrams with respect to both solid and waved lines (i. e., do not contain 
parts — ^ — and/or ^jxt^^jx, ). At large m most of the coupled skeleton 
diagrams of the (m+ 1)— th order can be constructed by coupling vertex >*\/w<; 
to a diagram ^J<HC^ of the m-th order with four outer lines obtained by breaking 
two of 2m solid lines in a coupled skeleton diagram of the m-th order. Since at 
large m the lines in most cases of coupled diagrams are arranged randomly (i. e., 
the symmetry elements are not important) this yields approximately 6m 2 (relevant 
to m ■ (2m — 1) possible ways to break the solid lines, and 3 different ways to couple 
the vertex to get a skeleton diagram) topologically nonequivalent diagrams of the 
(m + 1)— th order corresponding to a given diagram of the m-th order. However, 
different diagrams of the m-th order yield the same diagram of the (m + l)-th order. 
Really, we can subtract any of the m + 1 vertices ^>xna\><^ from a (random) 

coupled diagram of the (m + 1)— th order and obtain a diagram ^J^C^ of the m-th 
order, m + 1 reverse procedures represent m + 1 different ways of how to obtain this 
coupled diagram of the (m + l)-th order from diagrams of the m-th order. Thus, the 
total number of the above defined coupled skeleton diagrams of the (m+l)-th order 
exceeds the number of corresponding diagrams of the m-th order approximately 6m 
times, or the number of diagrams at large m increases approximately as 6 m m!. 

Let us consider other possibilities to construct coupled skeleton diagrams of the 
(m + l)-th order. We can couple a vertex ><\y\nj<^ to a diagram obtained by 
breaking two solid lines in a non skeleton diagram. The simplest of such diagrams 
^J^^ICi an d ^ifX^Ci are constructed of two above defined skeleton diagrams 
of orders m\ and m% with m\ + = m. The number of such diagrams relative to 
the number of the diagrams ^J^C^ obtained from the skeleton diagrams of the m- 
th order tends to zero at m — > oo as (wij mi!m2!) / (m 2 m!). More complicated 
diagrams like ^^Z^Z ••• Z^C^ (where any of the neighbouring blocks can be 
connected by waved lines instead of the solid lines) also can be included. However, 
this gives only a second-order correction. Other kind of diagrams among those ob- 
tained by breaking two solid lines in a non skeleton diagram, e. g., ^l^I^C^ , are 
not included because the coupling of vertex in this case does not yield a skeleton 
diagram. In such a way, the number of possible constructions is increased insignif- 
icantly if non skeleton diagrams of the m-th order are included. So, factor 6m is 
corrected only slightly. For a more accurate estimation this factor is replaced by 
6m + o(l), which allows to conclude that the number of diagrams of the m-th order 
in Eq. ( f29| ) increases with m like const ■ m 7 6 m m!, where 7 is a constant. Note 
that the diagrams of (|2^) are obtained from the actually discussed coupled skeleton 
diagrams by breaking a waved line, replacing this line by a (broken) dashed line. 

Consider now a random high-order diagram £ TO> j(q,£) of Eq. (p9|), i. e., the 
z-th diagram among those of the m-th order. Calculation of the sum over wave 
vectors and order-parameter-components for such a diagram is analogous to finding 
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of the statistical sum of a random lattice formed by the lines of this diagram. Two 
conditions must be fulfilled: first, the sum of the wave vectors coming into any 
of the kinks is zero, the wave vectors of the outer lines being fixed (q and — q), 
and, second, the same order-parameter-component is related to solid lines of a 
closed loop. These conditions represent certain interaction between the lines of the 
diagram. The "interaction energy" is zero if they are satisfied, and oo otherwise. 
Besides, each line has his ownw "energy". If normalized to T, this energy is equal to 
minus logarithm of the absolute value of the factor related to the corresponding line 
in the diagram notation. This analogy can be used in calculation of | £ m j(q, £) |. 
In a typical case of large m this yields S m) j(q, £) ~ — (— l) m b m , where b m with 
b > represents the absolute value and — (— l) m - the sign of this term. This result 
remains true (with new value of b) if the combinatorial factor is included, since the 



latter is 4-8 m for any random diagram having no symmetry elements (see Sect. 4.2 ) 



Appendix B. Non— perturbative analysis 

For convenience, here we consider the one-component case n = 1, since the 
specific (natural) value of n is irrelevant in the present scaling analysis. The Hamil- 
tonian of the considered model 

H/T = J [r </? 2 (x) + c(V9?(x)) 2 + up 4 (x)] dx (Bl) 

in the Fourier representation is 

H/T = ( r + ck 2 ) | ip k | 2 + uV' 1 WkiWsaPloV-ki-ka-ks ■ ( B2 ) 

k ki,k 2 ,k 3 

Summation in Eq. (B2) is performed within k < ko (here fco = A). Points in the 
k-space are separated by a distance 2tt/L, where L is the linear size of the system 
with volume V = L d at d < 4. By a formal change of variables = u a c~ da l 2 <{)\ L , 
p = c 2a u~ a h, po = c 2a u~ a ko, and R = roc da u~ 2a , where a = 1/(4 — d) we obtain 

H/T = J2(R + p 2 )\* P \ 2 + vr 1 £ ^ pi ^ P2 ^ p3 vi/_ pi _ P2 _ P3 (jB3 ) 

P P1.P2,P3 

where V\ = Lf with L\ = L-c~ 2a u a . In Eq. (B3) the summation is performed within 
p < Pq. Points are separated by a distance 2ir/Li . According to equations (B2) 
and (B3), we have (fluctuation modes with k > fco are excluded) 

2\ dct„—2a \ ,t, |2\ da„—2a 



G(k) = (| v?k r) = c aa u- 2a (| * p |^) = c aa u- 2a g(p,p ,R) (B4) 

where, at a fixed d and V\ — > oo, g is a function of the given arguments only. The 
latter is true since the only parameters contained in Hamiltonian (B3) are po, R, 
and V\. For the same reason, the critical value of R depends merely on po, and the 
correlation function at the critical point is 

G*(k) = c da u- 2a g*(p,p ). (B5) 
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It follows from Eq. (B5) that the width of the critical region for p is p cr u = (ftiPo) 
where <j>(z) is a single-argument function. Here both d and the accuracy parameter 
e in Eq. (^) are considered as fixed quantities. According to the relation between 
p and k, we have p cr it = c 2a u~ a k cr it, which yields 

kcrit = c~ 2a u a <P (c? a u- a k Q ) . (B6) 

Similarly, by using Eq. (B5) we obtain an exact scaling relation for the asymptotic 
expansion at k — ► 0, i. e., 

G*(k) = J2kk- X > , (B7) 
i 

where Xq = X = 2 — rj and terms with I > are corrections to scaling with the 
amplitudes represented in the scaled form 

b X = B t (c 2Q U- Q fc ) • C «( d - 2A ') n «(A i -2) _ (£g) 

It is evident from Eq. (B6) that k cr n can exponentially (like exp (— u~ a ) with a > 0) 
tend to zero at u — > only if the function 4>(z) decreases exponentially at z — > oo. 
Let us now consider the behavior of k cr a at k$ — > oo at fixed c and u. We conclude 
immediately: if k cr u decreases exponentially at u — * 0, then this quantity behaves 
in the same way at k$ — > oo. Physically, this means that the width of the critical 
region is strongly affected by short-wave fluctuations and the effect dramatically 
(exponentially) increases with shortening of the wavelength. The latter is rather 
unphysical, since the critical behavior is well known to be governed by long-wave 
fluctuations. Thus, k cr n cannot decrease exponentially (or, in general, faster than 

u s at any s > 0) at u — > 0, which means that lim (u s / k cr it(u)) = holds at large 

u— >o 

enough s. Based on Eq. (B8), the same conclusion can be made for the amplitudes 
k. 
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